Low-energy structures of zinc borohydride Zn(BH4)2 
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We present a systematic study of the low-energy structures of zinc borohydride, a crystalline ma- 
terial proposed for the hydrogen storage purpose. In addition to the previously proposed structures, 
many new low-energy structures of zinc borohydride are found by utilizing the minima-hopping 
method. We identify a new dynamically stable structure which belongs to the 74i22 space group as 
the lowest-energy phase of zinc borohydride at low temperatures. A low transition barrier between 
/4i22 and PI, the two lowest-lying phases of zinc borohydride is predicted, implying that a coexis- 
tence of low-energy phases of zinc borohydride is possible at ambient conditions. An analysis based 
on the simulated X-ray diffraction pattern reveals that the 74i22 structure exhibits the same major 
features as the experimentally synthesized zinc borohydride samples. 

PACS numbers: 61.66.-f, 63.20.dk, 61.05.cp 



I. INTRODUCTION 

Hydrogen is an environment-friendly cincrgy c;arricr 
which can provide a high energy density without pro- 
ducing greenhouse gases. This clean fuel is promising 
for many applications in, for example, the transporta- 
tion sector. Various modern generations of fuel-cell and 
hybrid vehicles are being commercialized, using hydro- 
gen as a fuel. Currently, hydrogen is stored in complex- 
structured high-pressure tanks of which the volume is a 
major difficulty for compact vehicles. A more efficient 
method for hydrogen storage is desired, and tremendous 
research efforts have been given to this goal^^^. 

Broad interest in metal borohydrides, a class of ionic 
crystal materials, is motivated by the possibility of using 
them for the hydrogen storage purpose^'^. Among these 
materials, alkali metal borohydrides are generally ther- 
modynamically too stable. For example, lithium boro- 
hydride LiBH4 can reversibly store 8 — 10% hydrogen at 
temperatures of 315 — 400°C, which is too high for on- 
board applications^. Zinc borohydride Zn(BH4)2, one of 
the divalent metal borohydrides, is an alternative because 
of its more favorable thermodynamical properties. In 
particular, Zn(BH4)2 has a low decomposition tempera- 
ture 85°C)''"^ and a relatively high gravimetric hydro- 
gen density (~ 8.5 wt %)^''^. However, the reversibility of 
Zn(BH4)2 remains poorly understood. While Zn(BH4)2 
was reported* to be reversible, an attempt to reduce the 
decomposition temperature and to enhance the kinetic 
by doping Zn(BH4)2 with Ni nanoparticles suppressed 
the reversibility^'^. Considerable research interest was 
therefore devoted to this crystalline material^""'^"^'"^^'"'^''. 

Experimentally, Zn(BH4)2 can be synthesized in sev- 
eral ways^"'^'^, one of which is through the metathesis 
reaction 

2NaBH4 + ZnCl2 -i> Zn(BH4)2 + 2NaCl. (1) 

Two X-ray powder diffraction (XRD) analysis for a 
Zn(BH4)2/NaCl mixture, the product of reaction (1), 



were reported in Refs. 6 and 7. However, the crystal 
structure of the Zn(BH4)2 products has not yet been con- 
clusively determined because of the insufficient chemical 
purity of the samples^'*". 

Detailed knowledge of the crystal structure is however 
essential for further studies of a variety of material prop- 
erties. Recently, ab-initio crystal structure prediction 
has become an increasingly attractive approach in ma- 
terial science, but still remains a challenging task (see 
Ref. [12] and references therein). Furthermore, density 
fimctional theory (DFT)^'^'^'' calculations have been ex- 
tensively used in the past to study hydrogen storage ma- 
terials, and its successes and shortcomings have been dis- 
cussed in detail by Herbst et al}^ . Several candidates for 
the low-temperature crystal structure of Zn(BH4)2 have 
been theoretically proposed in literature^^'^^'^^. All of 
the structures were predicted by the database searching 
method, starting from the existing crystal structures pro- 
posed for magnesium borohydride Mg(BH4)2, also a di- 
valent borohydride. A triclinic PI structure was initially 
suggested from the monoclinic P2/c structure predicted 
for Mg(BH4)2 in the same work by Nakamori et al}^ 
An orthorhombic Pmc2\ structure was subsequently pro- 
posed and shown to be dynamically stable by Choudhury 
et al}^ , starting from the Pmc2i structure initially in- 
tended for Mg(BH4)2^*. Two other structures, also orig- 
inally proposed for Mg(BH4)2, i.e., the tetragonal /4m2 
structure-'^^ and the orthorhombic F222 structure^^, were 
then examined by Aidhy and Wolverton in Ref. 16, 
yielding two low-energy nearly-degenerated structures for 
Zn(BH4)2. While both of the structures are much lower 
in energy (ci: 30 kJ mol~^ f.u."^) than the Pmc2i struc- 
ture, the F222 structure is energetically slightly favored 
over the /4m2 structure by ~ 1.0 kJ mol~^ f.u.~^ (here 
f.u. is used for "formula unit"). The F222 structure 
has therefore been used for the low-temperature phase 
of Zn(BH4)2 in a phase- stability analysis of some mixed- 
metal borohydride systems^^. 

In this paper, we revisit the low-temperature crys- 
tal structure of Zn(BH4)2 by first-principles calculations 
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TABLE I. Summary of the known structures (top panel) and the new structures (bottom panel) of Zn(BH4)2. Enthalpies of 
formation at OK {AHok) and lOOK (AT^iook) as well as their components {AHei, AHzp, and 5Ai?iooK) were obtained with the 
PRE functional. Total electronic energy ^p^EsoI ^^j^^j^out zero-point energy correction), obtained with the PBEsol functional, 

is given with respect to that of the J4i22 structure. The unit of the enthalpies and energies is kj mol~^ f.u.~^. Space group 
(SG) numbers are given within the parentheses nest to the corresponding space group symbols. 



Initial material Ref. Initial SG Final SG 


A77el 


A77zP 


<5A77iooK 


A77oK 


A77iooK 


^PBEsol 


Zn(BH4)2 


[16] 


F222 (22) 74i22 (98) 


-37.89 


3.11 


0.24 


-34.78 


-34.55 


0.00 


Mg(BH4)2 


[44] 


74i22 (98) /4i22 (98) 


-37.89 


3.11 


0.24 


-34.78 


-34.55 


0.00 


Zn(BH4)2 


[16] 


74m2 (119) 74m2 (119) 


-36.32 


- 


- 


- 


- 


1.64 


Zn(BH4)2 


[10] 


PI (2) PI (2) 


-26.45 


0.68 


0.06 


-25.77 


-25.71 


7.74 


Zn(BH4)2 


[11] 


Pmc2i (26) Pmc2i (26) 


-19.06 


3.41 


0.71 


-15.66 


-14.95 


18.50 






74i22 (98) 


-37.89 


3.11 


0.24 


-34.78 


-34.55 


0.00 






74m2 (119) 


-36.32 










1.64 






PI (1) 


-32.96 


2.23 


0.44 


-30.73 


-30.28 


5.92 






C2 (5) 


-32.71 


1.80 


0.29 


-30.91 


-30.62 


6.11 






Ibam (72) 


-30.47 


1.90 


0.29 


-28.57 


-28.28 


8.30 






C222i (20) 


-28.59 


-0.30 


-0.46 


-28.90 


-29.36 


14.92 






Ama2 (40) 


-27.44 


-1.49 


-0.86 


-28.92 


-29.78 


16.33 






Pm (6) 


-27.41 


-1.47 


-0.83 


-28.88 


-29.72 


16.18 



based on DFT. While all of the structures mentioned 
above are reexamined, we discover in addition a large 
number of low-energy structures of Zn(BH4)2 using the 
minima-hopping mcthod^^'^^. The transition barrier be- 
tween the two lowest-lying structures of Zn(BH4)2 is pre- 
dicted by means of concerted nudged elastic band calcula- 
tions. We then examine the dynamical stability of the ob- 
tained structures by phonon frequency calculations and 
analyze the simulated XRD patterns by comparing them 
with existing experimental results®'^. Finally, we discuss 
some relevant features of the crystal structure prediction 
methods used for Zn(BH4)2. 



II. COMPUTATIONAL METHODS 

First-principles calculations in this work were per- 
formed within the projector augmented wave formal- 
ism as implemented in the Vienna Ah Initio Simulation 
Package (VASP)^^"^^. We used the generalized gradient 
approximation with the Perdew-Burke-Ernzerhof (PBE) 
functionaP^ for the exchange and correlation energy. The 
semicore pseudopotential for zinc was used, of which the 
valence configuration is id^^As^. For boron and hydro- 
gen, the valence electron configurations are 2s^2p^ and 
Is^, respectively. The convergence of the total energy 
calculations was ensured by a 9 x 9 x 9 Monkhorst-Pack 
k-point mesh^'^ for sampling the Brillouin zone and a ki- 
netic energy plane wave cutoff of 800 eV. Atomic and 
cell variables were simultaneously relaxed until all the 
residual force and stress components were smaller than 
1 meV/A and 10~^ kbar, respectively. The space groups 
corresponding to the relaxed structures were determined 
by FINDSYM^*. 



For searching the low-energy structures of Zn(BH4)2, 
we used the minima-hopping method^^'^^, an efficient 
structure prediction approach which was recently ex- 
tended for crystalline systems. Different from other re- 
cent approaches which use, for example, the electrostatic 
energy ^^'^^ or the number of metal- hydrogen bonds^°, 
the minima-hopping method uses the energies evaluated 
at the DFT level as the objective function. The en- 
ergy landscape is explored by short consecutive molecular 
dynamics trajectories followed by local geometry relax- 
ations. The initial velocities for the molecular dynam- 
ics runs are chosen approximately along soft mode direc- 
tions, allowing efficient escapes from local minima, and 
aiming towards the global minimum. This method was 
successfully applied in a wide range of material structure 
predictions^^"^^, and is not at all restricted to ionic hy- 
dride materials but can be applied to any system of given 
compositions with arbitrary boundary conditions. Some 
of the theoretically predicted structures, e.g., the neu- 
tral Si clusters with more than 12 atoms and four- fold 
coordinated defect in silicon, were recently confirmed by 
experiments^^''*". 

The energetic ordering of the low-energy structures of 
Zn(BH4)2 was examined via the enthalpy of formation 
A77t, given by^^.^i 

AHt = AH^i + AHzp + SAHt- (2) 

Here, AHgi and AHzp are the electronic and zero-point 
energy differences between the products and the reac- 
tants, according to the reaction 1, while SAHt is the 
energy change from OK to T. Within the harmonic ap- 
proximation, AHzp and SAHt were straightforwardly 
computed from the results of the frozen-phonon calcula- 
tions as described in Section IIIB. We used the P4:2/mnc 



3 



phase of NaBH4''^, the Fm3m phase of NaCl, and the 
Pna2i phase of S—ZnCh'^^ for the calculations of AHt 
according to Eq. 1. 



III. LOW-TEMPERATURE STRUCTURE OF 
ZINC BOROHYDRIDE 

A. Low-energy structures 

We re-examined the crystal structures previously pro- 
posed by Refs. 10, 11, 16 for Zn(BH4)2 or by Ref. 44 
for Mg(BH4)2. Although the forces exerting on the ions 
of the F222 structure taken directly from Ref. 16 are 
quite smaU (< 0.05 eV/A), one still can gain ~ 1.7 kJ 
mol~^ f.u.~^ by further relaxation. The relaxed struc- 
ture was identified as an /4i22 structure, which can also 
be obtained by relaxing the /4i22 structure reported in 
Ref. 44 for Mg(BH4)2 after replacing the Mg atoms by 
Zn atoms. For the other structures, the corresponding 
space group remains unchanged after the relaxation. 

Several minima-hopping simulations were performed to 
search for additional low-energy structures of Zn(BH4)2. 
We found that the minima-hopping method is able to 
predict not only the existing /4i22 and /4m2 structures, 
but also a large number of new low-energy structures 
for Zn(BH4)2. We consider in this work the structures 
which were discovered within the energy range of ~ 10 
kJ mol~^ f.u.~^ above the lowest-energy structure. De- 
tailed information on these structures can be found in 
the supplemental material**^. 

The enthalpies of formation of all the examined struc- 
tures are shown in Table I together with their compo- 
nents, according to Eq. 2. At both OK and lOOK, 
the /4i22 structure (see Table II for detailed informa- 
tion) is the most thermodynamically stable structure of 
Zn(BIl4)2 . Table I also shows that the electronic en- 
ergy AHc\ of the /4to2 structure is slightly higher than 
that of the /4i22 structure by ~ 1.7 kJ mol"^ f.u."^ 
These two phases are energetically favored over the PI 
and Pmc2i phases by ~ 10 and by ~ 20 kJ mol~^ f.u.~^, 
respectively. These results, obtained with the PBE func- 
tional, are consistent with those obtained with the PW91 
functional for the exchange-correlation energy^^. 

To confirm the energetic ordering, we performed 
additional calculations for the electronic energy by 



TABLE II. Atomic positions of the tetragonal 74i22 structure 
for Zn(BH4)2. Cell parameters are a — b = 6.986A, c = 
12.189A, a = /3 = 7 = 90°. 



ABINIT"*^'^^. Our calculations were carried out 
with the norm-conserving Hartwigsen-Goedecker-Hutter 
pseudopotential'^®, a plane- wave cutoff energy of 60 Ha- 
tree (« 1600 eV), and the PBEsol functional, which is 
a modified PBE generalized gradient approximation that 
improves equilibrium properties of solids'*^ . The obtained 
results for the electronic energy £;PBEso1^ which are also 
shown in Table I with respect to that of the /4i22 struc- 
ture, is consistent with the energetic ordering obtained 
with the PBE exchange-correlation functional. 

The tetragonal /4i22 structure of Zn(BH4)2 is illus- 
trated in Fig. 1 (prepared by VESTA^°). Similar to the 
geometry of most of the other complex metal borohy- 
drides, the complex [BIl4]~ anions in the /4i22 structure 
for Zn(BH4)2 form isolated, slightly deformed tetrahedra 
with the B-H bond length of either 1.2lA or 1.24A while 
the H-B-H angle is either 104° or 119°. Each zinc atom 
is surrounded by four boron atoms with equal Zn-B bond 
length of 2.35 A and B-Zn-B angles of either 99° or 134° . 
The geometries of the other structures shown in Table I 
are somewhat similar, i.e., these ionic crystal structures 
are characterized by different arrangements of the Zn^+ 
cations and the complex [BIl4]^ anions in their more or 
less deformed tetrahedra. 



B. Structural stability 

The dynamical stability of the zinc borohydride struc- 
tures shown in Table I was determined by calculations 
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FIG. 1. (Color online) The tetragonal 74i22 structure for 
Zn(BH4)2. Cray (large), green (medium), and pink (small) 
spheres represent zinc, boron, and hydrogen atoms, respec- 
tively. 
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FIG. 2. (Color online) Calculated density of plionon states 
of the /4i22 and 74m2 structures for Zn(BH4)2. Imaginary 
plionon frequencies are represented on the figure by negative 
real values. 



of phonon frequencies using PHONOPY'^^ , a package based 
on the super-cell approach^^. For each relaxed structure, 
finite atomic displacements with an amplitude of 0.01 
A were introduced to a 2 x 2 x 2 super cell which con- 
tains 16 formula units of zinc borohydride (176 atoms). 
Calculations for the atomic forces within the super cells 
were then carried out by VASP, allowing for the second- 
order force constants to be determined^'^. The phonon 
frequencies of the structures were finally calculated from 
the dynamical matrices, given in terms of the force 
constants. The longitudinal optical/transverse optical 
(LO/TO) splitting was not taken into account since the 
effects of the LO/TO splitting were reported to be neg- 
ligible for other hydrides'^®'^^'^^. 

The obtained density of phonon states of the /4i22 
structure, which is shown in Fig. 2, implies that this 
structure is dynamically stable. We also found that the 
/4to2 structure is dynamically unstable because of two 
phonon modes, one of which presents at each of sev- 
eral high-symmetry points (see Fig. 2 for the density 
of phonon states). In particular, at Z, F, X, and N, the 
corresponding imaginary mode has a frequency of 2.41i 
THz, 1.79i THz, l.SOi THz, and 1.19i THz, respectively 
To explore these modes, we followed the corresponding 
atomic eigendisplacements at these points by the same 
procedure described in Refs. 55 and 56. The aforemen- 
tioned /4i22 structure was finally re-obtained by explor- 
ing the soft modes at either Z, TV, or X point. On the 
other hand, a dynamically stable structure which belongs 
to the 74 space group (# 82) was obtained by exploring 
the imaginary mode at F. However, the XRD pattern of 



the 14 structure, which is slightly lower than the /4m2 
structure by 0.2 kJ mol~^ f.u.^^, is identical with that 
of the IAm2 structure, clearly indicating that the Zn/B 
frames of these two structure are identical. Finally, the 
soft mode exploration of the IAm2 structure ended up 
with the dynamically stable /4i22 structure. 

The densities of phonon states of the other structures 
for Zn(BH4)2 are shown in the supplement material*^, 
indicating that they are dynamically stable. We note 
that the obtained density of phonon states for the Pmc2i 
structure, as shown in the supplement material^^, is con- 
sistent with the results reported in Ref. 11, which also 
indicates that the Pmc2i structure is dynamically stable. 



C. Structural transformation 

Polymorphism is commonly observed in molecular 
crystals and is driven by thermodynamically compet- 
ing phases separated by low transition barriers. Upon 
a closer examination of Table I we observe that a large 
number of structures exist in a small energy range of 
less than 10 kJ mol~^ f.u.^^ above the ground state. 
Therefore, coexistence or transformation of several dif- 
ferent phases at finite temperature might be possible, 
especially if phase transitions require only small activa- 
tion energies. A detailed and quantitative description 
of such a behaviour requires an accurate description of 
the free energy landscape, the kinetics of the transfor- 
mation, the melting temperatures of each phase under 
consideration, etc. Here we shall limit ourselves to es- 
timate an upper limit of the transition barrier between 
two selected dynamically stable phases, the triclinic PI 
phase and the tetragonal J4i22 phase, at K by em- 
ploying the generalized solid-state nudged elastic band 
(G-SSNEB) method" as implemented in the VASP TST 
tools. An initial concerted pathway was estimated by 
carefully selecting an appropriate representation of the 
unit cells, followed by an approximative nudged elastic 
band (NEB) simulation. The saddle point was then re- 
fined by employing the climbing image NEB (CI-NEB) 
approach until the gradients were converged to less than 
3 meV/A at the first order transition state. The ener- 
gies of the two end points and the saddle point structure 
were recomputed with higher accuracy using the HSE06 
hybrid functional^*"^° since PBE calculations are known 
to underestimates the barrier height, especially if the co- 
ordination number is reduced at the transition state®^. 

The CI-NEB transition pathway is shown in Fig. 3 to- 
gether with the structural evolution along the path. The 
planar structure in the PI phase is distorted and recom- 
bined to form a hexagonal network in the /4i22 phase. 
The barrier height was found to be 11.16 kJ mol~^ f.u.~^ 
(reactant) and 16.10 kJ mol^^ f.u.^^ (product), while 
with HSE06 we found 13.40 kJ mor^ f.u.^^ (reactant) 
and 19.47 kJ mol^^ f.u.~^ (product), respectively. These 
activation energies are roughly 2 — 3 times larger than 
the energy differences between the two phases and are 
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Reaction coordinate (a.u.) 



FIG. 3. (Color online) Left panel: structural transition from the triclinic PI (1) to the tetragonal 74i22 structure (4) for 
Zn(BH4)2 with intermediate steps. Gray (large), green (medium), and pink (small) spheres denote zinc, boron, and hydrogen 
atoms, respectively. Right panel: interpolated energies along the NEB pathway. Solid blue circles denote PEE results, whereas 
green crosses denote HSE06 energies. Structure at (2) corresponds to the saddle point. 



sufficiently low such that they could be easily overcome 
at ambient conditions. The phonon dispersion was cal- 
culated at the highest saddle point to confirm that a sin- 
gle imaginary phonon mode is present (see supplemen- 
tal material for details). A second mode was identified 
with partially imaginary frequencies close to the T point. 
This indicates that the saddle point itself has very low 
curvatures with phonons contributing to a high rate con- 
stant according to transition state theory. Although we 
only investigated this particular phase transformation we 
can assume that similarly low barriers can be found for 
structural transitions between other phases since all such 
transformations do not require the breaking of strong co- 
valent bonds but can be obtained by rearranging weakly 
bonded molecular subunits, as commonly observed in 
other molecular crystals. 



D. Structural identification 

While /4i22 was theoretically suggested to be the 
lowest-energy structure of Zn(BH4)2 at low tempera- 
tures, a comparison with available experimental data is 
particularly useful. Having some experimental XRD in- 
formation for Zn(BH4)2 at hand^'^ we were able to per- 
form such a comparative study. For this purpose, an 
XRD analysis was performed for all the examined struc- 
tures using FULLPRDF package^^. To be consistent with 
the reported experiments®'^, the Cu Ka radiation (wave- 
length A = 1.54 A) was used for the XRD simulations. 
The simulated XRD patterns of these structures, given in 



the supplemental material^^, demonstrate that the /4i22 
structure is different from the other structures, specifi- 
cally the i^222 structure. 

Our analysis on the simulated XRD patterns indicates 
that among the examined structures, the XRD pattern 
of the /4i22 matches favorably with the currently avail- 
able XRD data®'*". In Fig. 4 we show the simulated 
XRD pattern of the /4i22 structure together with some 
information extracted from the experiments. As pre- 
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FIG. 4. (Color onhne) Calculated XRD pattern of the /4i22 
structure for Zn(BH4)2. Solid diamonds indicate the major 
peaks that were not claimed to correspond to NaCl in Ref. 
6. Solid (open) circles indicate the orientations of the ma- 
jor peaks that were (was not) claimed to correspond to the 
experimentally synthesized Zn(BH4)2 sample in Ref. 7. 
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TABLE III. Total electronic energy E^^^ and dynami- 
cal stability of the examined structures for Zn(BH4)2 and 
Mg(BH4)2. The energy is given in unit of kj mol~^ f.u.^^ 
with respect to that of the /4i22 structure. 



Structure 


Zn(BH4)2 


Mg(BH4)2 


^PBE 


Stability 


^PBE 


Stability 


74i22 


0.00 


stable 


0.00 


stable 


Mm2 


1.68 


unstable 


0.10 


unstable 


C222i 


9.30 


stable 


57.33 


unstable 


Ama2 


10.45 


stable 


33.97 


unstable 


PI 


11.57 


stable 


31.27 


stable 


Pmc2i 


18.84 


stable 


21.69 


stable 



viously mentioned, two measured XRD patterns of the 
Zn(BH4)2/NaCl mixtures were reported in Refs. 6 and 7. 
In Ref. 6, several dominant peaks corresponding to NaCl 
were identified while some of the remaining unlabeled 
peaks are assumed to correspond to Zn(BH4)2. More re- 
cently, several peaks from 25° to 30° and from 35° to 45° 
were explicitly suggested to correspond to Zn(BH4)2 by 
Ref. 7. In Fig. 4, solid diamonds and solid circles indi- 
cate the positions of the peaks that were implied by Ref. 
6 and 7 to correspond to Zn(BH4)2. Although the infor- 
mation extracted from the experimentally observed XRD 
pattern is not sufficient for a conclusive structure deter- 
mination, a qualitative discussion on the identification of 
the examined structures for Zn(BH4)2 is possible. 

Fig. 4 indicates that the XRD pattern calculated for 
the /4i22 structure matches quite well with the XRD in- 
formation extracted from Refs. 6 and 7. There is a ma- 
jor peak located at the orientation of ~ 23° (indicated 
by an open circle), which was not identified to belong to 
Zn(BIl4)2 in Ref. 7 but can be seen in Ref. 6. At this 
orientation of the experimental XRD pattern by Ref. 7, 
there is however an intense peak identified to correspond 
to the sample holder. A possible overlap of the character- 
istic peaks of Zn(BH4)2 and the sample holder may be the 
reason explaining why this major peak of the Zn(BH4)2 
crystal was not explicitly identified^. 

Fig. 4 demonstrates that the /4i22 structure shares 
major structural features with that of the experi- 
mentally synthesized Zn(BH4)2 samples^'^. Obviously, 
the experimentally measured XRD patterns for the 
Zn(BH4)2/NaCl mixtures can not provide sufficiently ac- 
curate information for determining the low-temperature 
crystal structure of Zn(BH4)2. Therefore, XRD data for 
purified crystalline Zn(BH4)2 would be desirable for a 
proper determination of the structure. 



IV. CRYSTAL STRUCTURE PREDICTION 
METHODS FOR ZINC BOROHYDRIDE 

Prior to our investigations all the proposals for ionic 
crystal structures of Zn(BH4)2 were based on those of 



Mg(BH4)2^"'"'^^. It was also found^*' that the energy 
ordering of the F222, J4m2, and Pmc2i structures for 
Zn(BH4)2 is similar to the corresponding energy ordering 
of Mg(BH4)2, which phenomenologically shows a struc- 
tural correspondence between Zn(BH4)2 and Mg(BIl4)2. 
This observation is also supported by the so-called Gold- 
schmidt's rules of substitution^^, according to which one 
can substitute Mg atoms in an ionic crystal structure of 
Mg(BH4)2 by Zn atoms without disrupting the structural 
stability due to the similar charges and ionic radii of the 
Zn^+ and the Mg^+ cations. 

To further explore this structural correspondence, we 
determined the energy ordering and the dynamical stabil- 
ity of an extended list of structures, including three more 
structures obtained in this work, i.e., the PI, C222i, and 
Ama2 structures. For each of the structures, which were 
already obtained for Zn(BH4)2, Zn atoms were substi- 
tuted by Mg atoms, then the cell and the atomic vari- 
ables were fully relaxed. As reported by Ref. 44 and 
similar to the case of Zn(BIl4)2, the /4i22 structure for 
Mg(BH4)2 was also obtained by relaxing the F222 struc- 
ture. Phonon frequency calculations were then carried 
out to determine the dynamical stability of the struc- 
tures. A summary of this investigation is given in Table 
III and Fig. 5 while the densities of phonon states of these 
structures for Mg(BH4)2 are given in the supplemental 
material**^ . 

Table III and Fig. 5 show that, in agreement with 
Ref. 16, the energy ordering of the /4i22, /4to2, and 
Pmc2i structures for Zn(BH4)2 is similar to that for 
Mg(BH4)2. However, the energy ordering of the ex- 
tended list of structures for Zn(BH4)2 is different from 
that for Mg(BH4)2. Regarding the dynamical stability, 
four structures (/4i22, /4m2, PI and Pm,c2i) exhibit 
the same behavior for Zn(BH4)2 and Mg(BH4)2 while 
the other two structures {Ama2 and C222i) do not. Ob- 
viously, the phcnomenological structural correspondence 
between Zn(BH4)2 and Mg(BH4)2 is weak, especially for 
the low-energy phases which are closed in energy. 

It can be observed from Table I and also in Refs. 10, 
11, and 16 that the symmetry is not likely to break if a 
local geometry relaxation is started from a structure ob- 
tained by substituting Mg by Zn in a Mg(BH4)2 phase. 
Consequently, in principle, it is hard to explore the new 
symmetries which have not been reported in the litera- 
ture for Zn(BH4)2 and related materials, e.g., Mg(BIl4)2. 
This issue is however solved by an unconstrained search- 
ing method, e.g., the minima-hoping method. As illus- 
trated in Table I, one can easily explore the new symme- 
tries which have not been reported. 

It is worth noting that, while several previous 
studies^^'^^ reported that the Pmc2i structure for 
Mg(BH4)2 is dynamically unstable, the Pmc2i structure 
examined in this work for Mg(BIl4)2 was determined to 
be dynamically stable. To clarify this discrepancy, we 
have performed additional calculations for the Pmc2\ 
structure that was explicitly reported in Ref. 20. We 
found that, in agreement with the previous studies^^'^°. 
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Energy (kj/mol) 




10.45 
9.30 
1.68 
0.00 



I4m2 




I4m2 





14,22 




14,22 





57.33 

33.97 

31.27 

21.69 

0.10 

0.00 



Zn(BH4)2 



Mg{BH4 



FIG. 5. Energy ordering of the extended list of structures 
examined on Table III for Zn(BIl4)2 and Mg(BH4)2. Each 
structure is represented by a box with the corresponding 
space group. Numbers are given next to the boxes, indicating 
the energy difference (in kj mol~^ f.u.~^) between the cor- 
responding structures and the most stable structure. Lines 
connect the same structure for the two materials, indicating 
the re-arrangement of the examined structures when Zn atoms 
are substituted by Mg atoms. 



this structure is indeed dynamically unstable. On the 
other hand, the Pmc2i structure examined in this work, 
taken from Ref. 11 for Zn(BH4)2, is slightly deformed 
from that reported in Ref. 20, and using it for Mg(BH4)2 
yields a dynamically stable structure. For more infor- 
mation, the densities of phonon states of these slightly 
deformed Pmc2i structures for Mg(BH4)2 are shown in 
the supplemental material**^. 



V. CONCLUSIONS 

In conclusion, we have carried out a systematic study 
of the low-energy structural phases of Zn(BH4)2. By 
using the minima-hopping method, we have discovered 
many new low-energy structures of Zn(BH4)2. The most 
stable structure is identified to belong to the /4i22 space 



group. Phonon calculations demonstrate that the /4i22 
structure is dynamically stable while the /4m2 struc- 
ture, which is slightly higher in energy than the /4i22 
structure, is dynamically unstable. By following the 
atomic eigendisplacements corresponding to the unsta- 
ble phonon modes, the aforementioned dynamically sta- 
ble /4i22 structure is finally re-obtained. An XRD anal- 
ysis implies that the 74i22 structure shares some major 
structural features with the structure of the Zn(BH4)2 
samples that were experimentally synthesized. 

Furthermore, because many thermodynamically com- 
peting low-energy structures were found at low-energy 
during our structural search, and a concerted phase tran- 
sition between the two lowest-lying structures, i.e., the 
PI and the /4i22 structures, was investigated. We found 
that barrier connecting these two structures is low, im- 
plying that Zn(BH4)2 might exhibit a polymorphic be- 
havior at ambient conditions. 

The structural similarity between Zn(BH4)2 and 
Mg(BH4)2, mentioned in Ref. 16, plays an important 
role when proposing low-energy structures of Zn(BH4)2. 
We have discussed the similarity and found that, al- 
though there is a certain structural correspondence be- 
tween Zn(BH4)2 and Mg(BH4)2, it is weak and strongly 
limits the exploration of new symmetries of Zn(BH4)2 at 
low energies. We show that the minima-hopping method, 
on the other hand, allows for an efficient and fully un- 
constrained structural search of Zn(BH4)2. 
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TABLE I: Calculated parameters for the structures discovered by the 
minima-hopping method for zinc borohydride. For each structure, coll 
parameters are given while for each atom, the WyckofF site and the 
coordinates (x, y, and z) are given. 



PI 


a = 


7.37A, b = 6.96A, c = 6.42A 


(1) 


a = 


71.47°, P = 109.31°, 7 = 118.14° 


Zn 


(la) 


(-0.2109, -0.3116, 0.3580) 


Zn 


(la) 


(0.2494, 0.4200, 0.3545) 


TD 

D 


(^iaj 


(^U.UiU4, U.Udoo, U.zoOoJ 


B 


(la) 


fO 4fi98 -0 3169 0886') 


B 


(la) 


(0.0262, -0.4655, 0.4222) 


B 


(la) 


(-0.4210, 0.4392, -0.3726) 


H 


(la) 


(0.1871, 0.1170, 0.4143) 


H 


(la) 


(0.3214, 0.4931, 0.0885) 


H 


(la) 


(0.0100, -0.0349, 0.1638) 


H 


(la) 


(0.4313, -0.1989, -0.0870) 


H 


(la) 


(0.4251, 0.4745, -0.3780) 


H 


(la) 


(-0.4436, 0.3551, 0.4745) 


H 


(la) 


(0.1569, -0.3587, 0.3003) 


H 


(la) 


(0.0229, 0.3583, -0.4576) 


H 


(la) 


(-0.2764, -0.3701, -0.3687) 


H 


(la) 


(-0.3873, 0.3204, -0.1956) 


H 


(la) 


(0.0707, -0.3406, -0.4539) 


H 


(la) 


(-0.1452, 0.4872, 0.2852) 


H 


(la) 


(-0.1116, -0.0283, 0.4167) 


H 


(la) 


(-0.4959, -0.2305, 0.2439) 


H 


(la) 


(-0.0408, 0.2180, 0.1643) 


H 


(la) 


(-0.3831, -0.3565, 0.0822) 


C2 


fl< ^ 


6.98A, b = 12.98A, c = 6.29A 


(5) 


a 


90°, /3 = 71.67°, 7 = 90° 


Zn 


' 


(0.0000, -0.2392, 0.0000) 


Zn 


(2a) 


(O.OOOoi 0.0300, 0.0000) 


B 


(4c) 


(-0.3533, 0.3978, -0.2684) 


B 


(4c) 


(-0.2703, 0.1453, 0.0684) 


H 


(4c) 


(-0.4676, 0.3226, -0.2657) 


H 


(4c) 


(-0.2312, 0.2218, -0.0570) 


H 


(4c) 


(0.2193, 0.3964, 0.4442) 


H 


(4c) 


(-0.4433, 0.1316, 0.1885) 


H 


(4c) 


(-0.2804, 0.4003, -0.1153) 


H 


(4c) 


(-0.2342, 0.0719, -0.0648) 


H 


(4c) 


(0.4686, 0.4711, 0.2733) 


H 


(4c) 


(-0.1660, 0.1565, 0.1903) 


Prn 




6.74A, b = 8.00A, c = 5.67A 


(6) 


a = 


90°, 13 = 67.21°, 7 = 90° 


Zn 


(lb) 


(0.3547, 0.5000, -0.2154) 


Zn 


(la) 


(0.1220, 0.0000, -0.1081) 


B 


(lb) 


(-0.3953, 0.5000, 0.4045) 


B 


(2c) 


(0.2376, -0.2502, -0.0458) 


B 


(la) 


(-0.1276, 0.0000, -0.2582) 


H 


(lb) 


(-0.3596, 0.5000, -0.3896) 


H 


(2c) 


(0.1042, -0.1862, 0.1477) 


to be continued .. 



TABLE I — continued from previous page 



H 


(2c) (-0.1996, 0.1301, -0.2977) 


H 


(2c) (0.1301, 0.3438, -0.1254) 


H 


(2c) (0.3698, -0.3149, 0.0262) 


TT 
M 


(ibj (U.4Ud2, U.oUUU, U.44ooJ 


U 
XI 


(9o\ ( ^99Q ^^^QQ 9000^ 
t 1 i-u.o^zy, -u.OL>c/y, \j.Atjij\j\ 


XJ 

xl 


yJ.yjiLL^ U.UUUU, -yj.^yjZii.) 


XT 
XX 


l^ZCJ l^U.o'iD*, -U.IODU, -U.ZZOUJ 


XX 
XX 




C222i 


a = 11.85A, & = 6.65A, c = 7.15A 


(20) 


a = a = 7 = 90° 


Zn 


(4a) (0.3879, 0.0000, 0.0000) 


B 


(4a) (0.2061, 0.0000, 0.0000) 


D 


(4b) (U.UUUU, U.457U, U.ioUUj 


U 

LL 


t^oc J l^u.zuo^, u.iuuo, -u.uu^o^ 


XJ 
XI 


^^oc^ \_-u.^oyu, u.io/y, -u.iiou^ 


XJ 
XX 


^^oc^ ^^U-O^ioo, u.^iyou, -Kj.ou^xj) 


XJ 
XX 


I^OC j I^U.UOOO, U.'iOlO, -U.Z'iOOj 


Ama2 


1 5 — ; 7 1 

a = 8.00A, h = 12.44A, c = 5.67A 


(40) 


a = ^ = 7 = 90° 


Zn 


(4b) (0.2500, 0.0582, -0.1618) 


B 


(4b) (0.2500, 0.1831, -0.4266) 


B 


(4a) (0.0000, 0.0000, -0.0462) 


TT 
M 


(4bj (U.zoUU, U.zUiU, -U.zU44j 


XI 

XX 


l^oC J l^-U.UD'iO, -U.UOD'i, U.UoDl^ 


XJ 
XX 


{^n\ { n ^ftn9 n 9snQ n nnin^ 
l^oc j i^-u.oouz, u.zouy, u.uuiu^ 


XJ 
XX 


ypC) i^u.^uoi, u.^i^ioy, u.ozooj 


XX 
XX 


f^/iK^ /"n 9p;nn n ns'^R n/i7Qi'i 
l^'iDj (^u.zouu, u.uooo, -u.^/yij 


Ibam 


a = 7.68A, h = 9.25A, c = 6.84A 




u; — p — ^ — yu 


Zn 


(4b) (0.5000, 0.0000, 0.2500) 


B 


(8j) (0.3069, 0.0847, 0.0000) 


H 


(8j) (0.1559, 0.1189, 0.0000) 


H 


(8j) (-0.1836, -0.4524, 0.0000) 


H 


(16k) (0.1200, -0.3525, 0.1406) 



TABLE IL Crystallographic information of the PI saddle-point struc- 
ture on the CI-NEB transition pathway of Zn(BH4)2. Cell parameters 
are given while for each atom, the Wyckoff site and the coordinates {x, 
y, and z) are given. 



PI 


a = 


7.08A, b = 6.98A, c = 6.64A 


(1) 


a = 


78.12°, p = 111.71°, 7 = 113.67° 


Zn 


(la) 


(0.0724, -0.0936, -0.2237) 


Zn 


(la) 


(-0.3993, 0.2378, -0.3221) 


B 


(la) 


(-0.1833, -0.4233, -0.3321) 


B 


(la) 


(0.2446, 0.1748, -0.4899) 


B 


(la) 


(-0.2020, 0.0344, -0.2072) 


B 


(la) 


(0.3805, -0.0528, 0.0573) 


H 


(la) 


(0.0036, -0.3896, -0.2052) 


H 


(la) 


(0.1189, -0.0119, -0.4889) 


H 


(la) 


(-0.2015, 0.4616, -0.4623) 


H 


(la) 


(0.1535, 0.2913, 0.3999) 


H 


(la) 


(0.2130, -0.0229, 0.0522) 


H 


(la) 


(0.3561, -0.1333, -0.1055) 


H 


(la) 


(-0.0499, 0.1133, -0.2841) 


H 


(la) 


(-0.2205, -0.1271, -0.0944) 


H 


(la) 


(-0.4794, 0.1212, 0.0677) 


H 


(la) 


(0.4087, -0.1771, 0.2109) 


H 


(la) 


(-0.1694, 0.1680, -0.0887) 


to be continued .. 
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FIG. 1: (Color online) Calculated density of phonon states of the PI, Pmc2i, PI, C2, Ibam, C222i, Ama2, and Pm structures 
for Zn(BH4)2. 
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(la) (-0.3689, -0.0050, -0.3690) 
(la) (-0.3007, -0.4975, -0.2140) 
(la) (0.3376, 0.2399, -0.3043) 
(la) (-0.2342, -0.2787, -0.4424) 
(la) (0.3669, 0.1636, 0.4207) 



4 




FIG. 2: (Color online) Calculated density of phonon states of the highest (PI) saddle point structure on the CI-NEB transition 
pathway of Zn(BH4)2. See Fig. 3 of the text for more information. 
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Frequency (THz) 

FIG. 3: (Color online) Calculated density of phonon states of the /4i22, J4m2, PI, Pmc2\, C222i, and Aina2 structures for 
Mg(BH4)2. 



5 




FIG. 4: (Color online) Density of phonon states of the Pmc2i structures for Mg(BH4)2 taken from Voss et al, J. Phys. Condens. 
Matter 21, 012203 (2009) — solid curve, and from Choudhury et al., Phys. Rev. B 77, 134302 (2008) — dashed curve. 
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FIG. 5: (Color online) Calculated XRD patterns of the dynamically stable structures for Zn(BH4) 



